Wake effects characterization using wake oscillator model 
Comparison on 2D response with experiments 



BenoTt Gaurier 1 *, Gregory Germain 1 and David Cebron 2 

1 Ifremer - Bassin d'Essais, 150 quai Gambetta - BP 699 
62 321 Boulogne- sur-Mer, France 

2 Ecole Centrale Nantes, 1 rue de la Noe - BP 92101 
44 321 Nantes Cedex 3, France 

* Corresponding author: benoit.gaurier@ifremer.fr 



A model using wake oscillators is developed to predict the 2D motion in a transverse 
plan of two rigid cylinders in tandem arrangement. This model of the wake dynamics 
is validated with experimental data from previous trials which took place at the 
Ifremer flume tank in Boulogne-sur-Mer, France. The agreement between the model 
and the experimental results allows using this model as a simple computational tool 
in the prediction of 2D Vortex-Induced Vibrations (VIV) and, after some futher 
developments, Wake-Induced Oscillations (WIO) effects. 

1. Introduction 

Mooring and flow lines involved in offshore systems for oil production are submitted to various 
solicitations. Among them the effects of current are dominating. Vortex- Induced Vibrations (VIV) 
and Wake-Induced Oscillations (WIO) on closely spaced marine risers may lead to fatigue, clashes 
and structural failures. Extended studies have been conduced to describe and explain them for 
spring mounted uniform cylinders in translation perpendicularly to their main axis [H]. In the case 
of a pivoted cylinder with uniform diameter [6] a similar response is observed with some variations 
on the reduced velocity interval and maximum response. 

Experiments on model scaled tests with real configurations for dual risers interaction in a uniform 
current were performed in the Ifremer flume tank, within the framework of the project Clarom 
CEPM CO 3007/04 in partnership with Doris engineering, Saipem S.A., Institut Frangais du 
Petrole, Oceanide, Ecole Centrale Marseille and Total. The behaviour of two risers exposed to 
steady current and activated by VIV and WIO were studied [8] . These tests give a lot of informa- 
tion on how fluid interaction between two cylinders of equal diameter in tandem configuration can 
significantly modify their structural response in term of amplitude and frequency, compared to a 
single one [7]. Both in-line and cross-flow response have been studied and presented as functions 
of the reduced velocity. Results demonstrate that wake effects can be relatively strong. In almost 
all the tested cases the upstream cylinder responds like an isolated single one, whereas the vortex 



shedding and synchronization of the downstream cylinder can be strongly affected by the wake of 
the upstream one. Those phenomenons are relatively hard to predict. 

In order to quantify those wake effects, we developed a 2D phenomenological model of the near 
wake based on Van Der Pol wake oscillator (see [5] and [16] ) which describes the 2D motion of the 
cylinder in its transverse plan. This simplified model of the wake dynamics was first validated on 
a single cylinder in [3J and is extended here for the case of two cylinders in interaction. 

2. Mathematical model 

2. A. Dynamic equations in 2D 

The equations of this mathematical model are presented in [JJ. Here are just the main formulas 
needed for the basic comprehension. 

Initially, dynamic equations are written in a generalised way for a cylinder in free motion in its 
transverse plan. For an oscillating cylinder (mass m, volume V, velocity X) submit to a current 
described by: U = U x (t) + U y (t) e^, the center of inertia theorem gives: 

ma}- = F<^t + fie + Tie (1) 

with a£ = the relative acceleration, /j e = — m U the inertial force and fi c = the Coriolis 

force. The exterior forces Fext contain the hydrodynamic forces (drag, lift and forces issued from 
the potential theory), the spring force and possibly structural damping forces: 



^ Fext — Fhydro ~i~ f spring ~i~ f damping (2) 

Those previous forces are those in laminar flow, without vortex. But in turbulent flow, other 
special forces have to be added: the fluctuating drag and lift forces created by vortex and the 
blockage drag which is the additional drag issued from the transverse motion of the cylinder: this 
motion increases the apparent projected surface in front of the flow. 

Finally, in projection with Fl = et x Fp, the final general equations are: 

(m + C mP V)x + Xx + k(x - x ) = -pS [{C D + C Df ){U x - x) - (C L + C Lf )(U y - y)} x 

1 



y/(U x - ±) 2 + (U y - y) 2 + pV(l + C m )U x + - P SC Db ^U 2 + U* U x 

(3) 

(m + C m pV)y + Xy + k(y - y ) = -pS [(C D + C Df )(U y - y) + (C L + C Lf )(U x - x)] x 

1 



y/(U x ~ X) 2 + (Uy ~ ij) 2 + P V(1 + C m )Uy + - P SC Dh ^JU 2 + U 2 Uy 

with p the mass density, C m the added mass coefficient, Co and Cpf the average and fluctuating 
drag coefficient, Cob the blockage drag coefficient, Cl and Cl/ the average and fluctuating lift 
coefficient, k the stiffness and A the linear structural damping. 



2.B. Vortex forces 

Following [5], the vortex forces could be modeled on Van Der Pol oscillators coupled with the 
acceleration of the cylinder: 



( 

( 




- lj C Df + (2oj Cl ) 2 C Df = A D (x - U x ) 
l)d Lf + u 2 CL C Lf = A L (y-Uy) 



(4) 



where Cjjf and Clj q are the amplitudes of the fluctuating drag and lift coefficients, oj st is the 
vortex shedding pulsation, issued from the Strouhal number. [5] and [16] use Al = 12 and el = 0.3 
for one cylinder in ID motion, while we use here: Al = 7, Ad = 1, ££> = 1.2 and el = 2.5 for 
both cylinders. These parameters are issued from a preliminary study based on an optimization 
algorithm: the method of gradient descent. 

2. C. Hydrodynamical parameters 

For the accuracy of the model, the evolution of hydrodynamical coefficients must be known precisely. 
So, bibliographical data from |13j . |14j . |15j . |11] and [3] is fitted to obtain analytic formulas, versus 
Reynolds number Re = UD/v or reduced velocity V r = U/(f n D), with f n the natural frequency of 
the system. All these formulas are presented in [1] and concern the mean drag coefficient Co, the 
Strouhal number St, the added mass coefficient C m , the correlation length Al and the fluctuating 
lift coefficient C^f- 

Indeed, these formulas refer to only one cylinder. For the downstream cylinder, there are some 
differences to take into account, linked to the presence of the upstream one and of its wake. For 
example, on the contrary to the upstream cylinder, the average lift force of the downstream cylinder 
Cl 2 is not always equal to zero when this cylinder is not in line with the upstream one. 
In the following lines, indice 1 stands for the upstream cylinder, whereas indice 2 stands for the 
downstream cylinder. 

For considering these wake effects, [2] proposes to use the formula of the velocity deficit, first 
established by [15], and adjusted to experiments (figure [T]): 



with a\ = 1.2 and 02 = 13. 

According to [2], the equation [5] is valid when x is more than a few cylinder diameters. The 
hypothesis is that the drag force on the downstream cylinder is reduced from that in the free 
stream because the velocity in the wake is lower than the free stream velocity. So, it seems to be 
reasonable to estimate this force on a cylinder in a wake with the local wake velocity: 
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Fig. 1. Velocity deficit behind a cylinder, according to 



CD 2 (x,y) = C Dl 



Cdi 



,C Dl D / y 
1 — a\\j exp -a 2 



C Dl D x 



(6) 



with, according to [2J, a\ = 1 and a 2 = 4.5. 



By admitting Price's hypothesis [12], which postulate that the lift is proportional to transverse 



gradient of drag Cl 2 cx 



dCp 2 
d(y/D) 



the lift coefficient could be given: 



Cl 2 {x,v) = a 3 



C Dl x 



Cd 2 {x,v) y 



, C Dl D 
1 — ai\j exp -a 2 



C Dl Dx 



C Dl D 



exp —a 2 



C Dl D x 



(7) 



with, according to [2J, a\ = 1, a 2 = 4.5, as before, and 03 = —10.6. 



An implicit hypothesis is made here: these formulas, established for fixed cylinders, are supposed 
to be true for each time step. So the drag and lift on the downstream cylinder have to be calculated 
for each time step, especially because of the time-dependent distances x = x 2 — x\ and y = y 2 — yi 
between the two cylinders. 

In addition, parameters like C m2 , 0J s t 2 , Cif 2 and Cob 2 have to be calculated with the local velocity 
from equation For the stability of the code, this velocity is considered constant and calculated 
from the initial position of the upstream cylinder. 

Moreover, experimental data has shown that from a certain velocity, the motion of the down- 
stream cylinder is chaotic and the oscillation frequency which is used in the calculus of the Cob 2 is 
not really defined. In fact, the transverse motion of this cylinder is not periodic enough for gener- 
ating an additional drag. This phenomenon is directly linked by the turbulence on the wake of the 
upstream cylinder. So, with these considerations, the blockage drag of the downstream cylinder is 
chosen equal to zero. 



2.D. Link with experimental data 

To compare this model with our experimental data [7j, we have to link it to a pendulum motion of a 
rigid cylinder elastically mounted in a flow. Considering little angles, we can transform the pendulum 
equations in linear translation equations. To do this, we have to use, in the linear equation, the 
mass and the stiffness linked with the slice of the cylinder considered and located at z = 1.348 m. 
We use: 

I K 

m = and k = (8) 

with / the moment of inertia, K the angular stiffness and m the equivalent mass and k the equivalent 
stiffness. We use a structural damping coefficient £ = 1% in agreement with experimental free decay 
test in calm water. 

2. E. Algorithm 

The numerical scheme used to solve this problem is the implicit numerical differentiation formulas 
ode 15s of orders 1 to 5, from Matlab and especially designed for stiff systems PQ. 

To solve equations [3] and HI we need also to know Cr>b and C m , which are dependant on the 
amplitude A and of the frequency of transverse oscillation f ex . To solve this dependency problem, 
the algorithm used makes iterative loops while Cob{t + alt) — Cob{t) and C m (t + dt) — C m (t) are 
greater than a certain value, fixed here at 0.01. Then, when C^b and C m are converged, the stream 
velocity is incremented and the motion achieved. 

3. Comparison with experiments 

First, a free decay test is done in order to check the natural frequency calculated and to give a 
initial validation of the model for one cylinder [3J. The result is not presented here, but shows an 
agreement for both frequency and amplitude between the model and the experimental results. 

The mean hydrodynamical coefficients calculated by the model during this test are: C m = 0.98, 
St = 0.215 and Crjb = 0- These values correspond to the classical bibliographical results. 

3. A. Free cylinder in a flow 

To validate the model in a large range of reduced velocity, the complementary characteristics: mean 
and standard deviation displacements can be compared between experiments and model results. 



On figure 2(a) the mean transverse displacement of the upstream cylinder is of course null for all 
the reduced velocities, contrary to the mean in-line displacement, which is always increasing with 
the velocities. The change on slope at V r = 10 comes from the end of the lock-in. This sudden gap 
appears at a lower reduced velocity for the model (at V r = 10) than for experiments (at V r = 11). 



On figure 2(b) , the three branches introduced by [9] are plotted on the transverse displacements: 
the initial branch for Vr < 6, the upper branch for 6 < V r < 10 and the lower branch for V r > 10. 
Comparing with the experimental data, the gap between the upper branch and the lower branch 
appears at V r = 10 for the model, whereas it appears at V r = 11 experimentally. This difference 



(a) mean (b) standard deviation 



Fig. 2. Displacements of the upstream cylinder. Transverse oscillations: A model, 
A experiments. In-line oscillations: ■ model, □ experiments. The dash line is the 
quasi-static result. 



comes from the instabilities observed during experiments for these velocities: the transition from 
the upper to lower branch involves an intermittent switching between two modes [7]. 

Finally, in-line r.m.s. displacements seem to be less in agreement with the experimental data than 
the previous data. However, the comparison with bibliographical results like [ID] shows that the 
model reproduce correctly the in-line r.m.s. displacements for V r < 4, including the amplification 
of the motion at 1.7 < V r < 2.3 named the second instability region by [14] , Model results are also 
in relative agreement for V r > 10. However, these displacements are not reproduced during the 
lock-in. The origin of this difference comes from the amplification of the amplitude of Crjf 01 which 
is not take into account here. Indeed, contrary to the amplitude of the fluctuating lift coefficient 
Clj 01 , no bibliographical data was found on this subject. So, this coefficient is considered constant 
in the model. 



For the case presented here, the downstream cylinder is initially placed at 5D in line behind the 
upstream cylinder. On figure 3(a) , the mean transverse displacement of the downstream cylinder 
is also always null because the two cylinders are initially in tandem arrangement. The mean in-line 
displacement is increasing monotonically with the velocity but shows a short plate around V r = 10 
when the upstream cylinder's lock-in stops. As explained in paragraph s. CI there is no amplification 
of the mean drag, and so of the mean in-line displacement: the slope corresponds to the quasi-static 
curve. The small differences observed for 12 < V r < 14 and V r > 18 come probably from the 
differences noticed in the in-line r.m.s. displacements of the upstream cylinder. 

On the figure |3(b)| the numerical results are quite different from the experimental ones. The 
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Fig. 3. Displacements of the downstream cylinder. Transverse oscillations: ▲ model, 
A experiments. In-line oscillations: ■ model, □ experiments. The dash line is the 
quasi-static result. 

standard deviation of the displacement is in agreement especially for 8 <V r < 14. For the highest 
velocities, the wake instabilities are rather well represented by the transverse motions, but not for 
the in-line displacements. This difference comes essentially by the lack of data on the fluctuating 
amplitude of the drag force Cjjf 01 and Cr)f 02 for both cylinders. It could be very interesting to make 
new experimental trials to measure these forces on two fixed and rigid cylinders, with many initial 
positions, in order to determine precisely the evolution of these forces in function of the relative 
position of the cylinders. 

4. Conclusion 

The behaviour of two cylinders in tandem arrangement at 5D and submit to a current had been 
modeled by a phenomenological model based on Van Der Pol oscillators. After the presentation of 
the model and the description of the used parameters, we have compared the model results with 
some experimental data. The comparison is relatively in agreement for the mean and standard 
deviation of the two dimensional cylinders motions. The differences come essentially from the poor 
set of data for the fluctuating drag coefficient. Some specific experimental tests should be conducted 
to solve this problem. 

Despite those imperfections, other initial configurations will be tested and compared between 
experimental data and this simple model, to extend and determine the limitations of this kind of 
two dimensional code. 
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